Decay of an excited atom near an absorbing microsphere 
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Spontaneous decay of an excited atom near a dispersing and absorbing microsphere of given com- 
plex permittivity that satisfies the Kramers-Kronig relations is studied, with special emphasis on a 
Drude-Lorentz permittivity. Both the whispering gallery field resonances below the band gap (for a 
dielectric sphere) and the surface-guided field resonances inside the gap (for a dielectric or a metal- 
lic sphere) are considered. Since the decay rate mimics the spectral density of the sphere-assisted 
ground-state fluctuation of the radiation field, the strengths and widths of the field resonances essen- 
tially determine the feasible enhancement of spontaneous decay. In particular, strong enhancement 
can be observed for transition frequencies within the interval in which the surface-guided field reso- 
nances strongly overlap. When material absorption becomes significant, then the highly structured 
emission pattern that can be observed when radiative losses dominate reduces to that of a strongly 
absorbing mirror. Accordingly, nonradiative decay becomes dominant. In particular, if the distance 
between the atom and the surface of the microsphere is small enough, the decay becomes purely 
nonradiative. 

PACS numbers: 42.50.Ct, 42.60. Da, 42.50.Lc 



I. INTRODUCTION 

Light propagating in dielectric spheres can be trapped 
by repeated total internal reflections. When the round- 
trip optical paths fit integer numbers of the wavelength, 
whispering gallery (WG) waves are formed, which com- 
bine extreme photonic confinement with very high qual- 
ity factors. The frequencies and linewidths of WG waves 
are highly sensitive to strain, temperature, and other pa- 
rameters of the surrounding environment in general. The 
unique properties of WG waves are crucial to cavity- 
QED and various optoelectronics applications WG 
waves with Q- values larger than 10 9 have been reported 
for fused-silica microspheres (2|-Q| and liquid hydrogen 
droplets Q , and the ultimate level determined by intrin- 
sic material absorption has been achieved ^ . 

Since the spontaneous decay of an excited atom de- 
pends on both the atom and the ambient medium ||, it 
can be expected that if an atom is near a dielectric mi- 
crosphere, its spontaneous decay sensitively responds to 
the sphere-assisted electromagnetic field structure. In the 
case of an atom in free space, the well-known continuum 
of free-space density of modes of radiation is available 
for an emitted photon. By redistributing the density of 
possible field excitations in the presence of dielectric bod- 
ies, the radiative properties of excited atoms can be con- 
trolled. Experimental observations of lifetime modifica- 
tions of ions or dye molecules embedded in microspheres 
or droplets 0, cavity QED effects in the coupling of a 
dilute cesium vapor to the external evanescent field of a 
WG mode in a fused silica microsphere || , and detection 
of individual spatially constrained and oriented molecules 
on the surfaces of glycerol microdroplets || have been 
reported. In Refs. ju^ll]], the decay rate of a two-level 
atom located in the vicinity of a dielectric sphere was cal- 
culated and enhancement factors of hundreds were pre- 



dicted. The calculations are based on the assumption 
that there is no material absorption. In practice how- 
ever, there exists always some material absorption that, 
even when it is small, can still affect the system in a 
substantial way. Knowledge of both the radiative losses 
and the losses due to material absorption is crucial, e.g., 
in designing ultralow threshold microsphere lasers using 
cooled atoms or single quantum dots |l2| . 

In the present paper we study the problem of sponta- 
neous decay of an excited atom near an absorbing mi- 
crosphere of given complex permittivity that satisfies the 
Kramers-Kronig relations, applying the formalism devel- 
oped in Refs. |D|,[l4|]. This enables us to take into ac- 
count both dispersion and absorption in a consistent way, 
avoiding any restrictive conditions with respect to the fre- 
quency domain. In the numerical calculations we assume 
that the permittivity of the microsphere depends on the 
frequency according to the Drude-Lorentz model, which 
covers dielectric and metallic matter. 

The paper is organized as follows. In Sec. [n] the radi- 
ation field resonances associated with a microsphere are 
examined. Equations for determining the positions and 
widths of both the WG field resonances below a band 
gap and the surface-guided (SG) field resonances inside 
a band gap are given. The corresponding quality factors 
are calculated and the effects of dispersion, absorption, 
and cavity size are studied. Further, the spatial distribu- 
tion of the cavity-assisted ground-state fluctuation of the 
radiation field is studied. In particular, the wings out- 
side the microsphere of the fluctuation of the WG and 
SG waves are just responsible for exciting such waves in 
the spontaneous decay of an excited atom that is situated 
near the surface of the sphere. 



In Sec. HI the rate of spontaneous decay, the line shift , 
and the (spatially resolved) intensity of the emitted light 
are calculated as functions of the atomic transition 
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FIG. 1. Real and imaginary parts of the permittivity of 
a Drude-Lorentz-type dielectric, Eq. (lj), for up — 0.5 lot and 
7 — 10" 4 lut- The band gap covers the interval from lot to 
ljl — 1.12 LOT- 

frequency, with special emphasis on the influence of the 
orientation of the transition dipole moment, the distance 
between the atom and the microsphere, and the disper- 
sion and absorption outside and inside a band gap. In 
this context, the ratio of radiative decay to nonradiative 
decay due to material absorption is analyzed. 

In Sec. |fy] the results derived for the general case of 
dielectric sphere are applied to a metallic microsphere by 
appropriately specifying the matter parameters. Typical 
features are briefly discussed. Finally, some concluding 
remarks are given in Sec. 0. 



II. RADIATION FIELD RESONANCES 

Let us consider a dielectric sphere of radius R whose 
complex permittivity e(ui) = en(ui) + iej(uj) reads, accord- 
ing to a single-resonance Drude-Lorentz model, 
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(1) 



where top corresponds to the coupling constant, and lot 
and 7 are respectively the medium oscillation frequency 
and the linewidth. An example of the dependence on 
frequency of the permittivity is shown in Fig. |l|. Note 
that e(w) satisfies the Kramers-Kronig relations. From 
the permittivity, the complex refractive index can be ob- 
tained according to the relations 



n{uo) = \J e(u) — u r (lo) + inj(uo), 
™i?(/)M = i 



V^)+¥) + (-)eflH 



(2) 
(3) 



The dielectric features a band gap between the trans- 
verse frequency lot and the longitudinal frequency 
Far from the medium resonances, we 
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typically observe that 



e/(w) < \e R (oj)\. (4) 
For uj < lot, i-e., outside the band gap, we have 

Cfl(w) > 1, (5) 

n R (oj) ~ ^Je R {u) > ni(uj) ~ ei }-l= > (6) 

2^e R (Lo) 

and for lot <lo < lol, i.e., inside the band gap, 

< 0, (7) 
u r (lo) ~ 6I{UJ) « ~ V\^rW\- (8) 

VM W )I 

Note that ci(lo) > in any case. 

Using the source-quantity representation of the quan- 
tized electromagnetic field Q,^5|, all the information 
about the influence of the sphere on the spontaneous de- 
cay of an atom is contained in the Green tensor of the 
classical, phenomenological Maxwell equations, where 
dispersion and absorption are fully taken into account 
In particular, the sphere-assisted (transverse) 
field resonances that can give rise to an enhancement of 
the spontaneous decay are determined by the poles of the 
transverse part of the Green tensor. 

From the explicit form of the Green tensor as given in 
Appe ndix A, i t foll ows, o n se tting the d enominators in 
Eqs. @ [or ([AH])] and (|a|) [or jAlT| )1 equal to zero, 
that the resonances are the complex roots 



to = Q — id 



of the characteristic equations 



M(lo) = 



H' l+1/2 (hR) 
H l+1/2 (hR) 

for TE waves, and 
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J' l+1/2 (k 2 R) 
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1 - 



= (11) 



for TM waves, where k\ — lo/c, k 2 = \J e(lo)lo/c, and R 
is the microsphere radius [Ju(z) - Bessel function; H u (z) 
- Hankel function]. In Eq. (g), Q and S, respectively, 
are the position and the HWHM of a resonance line. 
Equations ( |Io|) and ([ll]) arc in agreement with the re- 
sults of classical Lorenz-Mie scattering theory |l| . Similar 
equations with real permittivity have also been derived 
quantum mechanically by expanding the field operators 
in spherical wave functions |l(| . Note that these results 
cannot be extended to absorbing media by simply replac- 
ing the real permittivity by a complex one, but requires a 
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more refined approach to the electromagnetic fiel d q uan- 
tization in absorbing dielectrics (for details, see Jjjl and 
references therein). 

The linewidths of the resonances are determined by 
radiative losses associated with the input-output cou- 
pling and losses due to material absorption. For small 
linewidths, the total width can be regarded as being the 
sum of a purely radiative term and a purely absorptive 
term. In that case we may write 



Jrad 
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(12) 



where S rac i corresponds to the linewidth obtained when 
the imaginary part of the permittivity is set equal to zero. 



A. Resonances below the band gap 

We first consider the (transverse) field resonances be- 
low the band gap, ui < cut, where WG waves can be ob- 
served. These resonances are commonly classified by 
means of three numbers @^]: the angular momentum 
number I, the azimuthal number m, and the number i 
of radial maxima of the field inside the sphere. In the 
case of a uniform sphere, the WG waves are (21 + l)-fold 
degenerate, i.e., the 21 + 1 azimuthal resonances belong 
to the same frequency Q;^. 

When the imaginary part of the refractive index is 
much smaller than the real part, ni(uj)<^n R (uj), then the 
method used in Ref. JllJ for solving Eqs. ( |Io|) and (|ll|) in 
the case of constant, real refractive index formally applies 
also to the case of frequency-dependent, complex refrac- 
tive index. Interpreting the WG waves as resulting from 
total internal reflection, it follows that RRc fc 2 v > Rki , 
where v = l + 1/2, and v can be assumed to be large in 
general. Following Ref. [jl7|, the complex roots ui t i of 
Eqs. © and (0) can then formally given by 



= fi,i[n(u)i,i)] + 0(v 1 ), 



(13) 



where 
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with P = n(uiij) and P = l/n(u)i t i), respectively, for TE 
and TM waves, and the a, are the roots of the Airy func- 
tion Ai(— z). Although Eq. ( |l3|) is not yet an explicit 
expression for the roots as in the case of a nondispersing 
and nonabsorbing (idealized) medium, it is numerically 
more tractable than the original equations and offers a 
first insight into the WG waves. 



1. Frequency-independent permittivity 

In many cases in practice, it may be a good approxi- 
mation to ignore the dependence on frequency of the re- 
fractive index within a chosen frequency interval. From 
Eqs. ( |I3| ) and (|l4|), the positions of the WG waves are 
then found to be 



n = f(n R ) + o 



n R 



(15) 



For notational convenience, here and in the following we 
drop the classification indices. Obviously, the first term 
in Eq. ( [l5|) provides a very good approximation for the 
positions of the WG waves, provided that the ratio nijn R 
is sufficiently small. Note that n~1.45-M10 -11 for silica 
H and n ~ 1.47 + il0~ 7 for glycerol Q in the optical 
region. In other words, small material absorption does 
not affect much the positions of the WG waves. 

The contribution to the linewidth of the material ab- 



sorption can be evaluated from Eq. (13) by first order 
Taylor expansion of f(n) at n R , 



Sabs - -nif'(n R ) 



(16) 



To calculate the total line width determined by both 
radiative and absorption losses, one has to go back to 
the original equations (10) and (pi]), since the radia- 
tive losses are essentially determined by the (disregarded) 
term 0(i/ _1 ) in Eq. (p"3|). First-order Taylor expansion 
of M(w) around f2 yields 



Im 



M(Q) 



M'(fi) 



where 



M'(Q) 



for TE waves, and 

M'(fi) ~ (e - 1) 
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(17) 



(18) 



(19) 



for TM waves. The radiative contribution to the line- 
width can then calculated by means of Eq. (|l2|). 

From Fig. || it is seen that the contribution to the 
linewidth of the material absorption increase with the 
imaginary part of the refractive index nj and the radius 
R of the sphere. It is seen that when the values of ri] and 
R rise above some threshold values, then the absorption 
losses start to dominate the radiative losses. It is fur- 
ther seen that for chosen nj the threshold value of R 
decreases with increasing n R . For instance, the absorp- 
tion losses start to dominate the radiative losses if R > 8A 

j(a)], and if R>1.6\ for 
It should be noted, that 



for n= 1.45 + HQ 9 [see Fig. 
n = 2A5 + ilQ- 9 [see Fig. |(b) 



the approximate results based on Eqs. ( |l6| ) and (|l7|) are 
in good agreement with the exact ones (without disper- 
sion) . 
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2. Frequency-dependent permittivity 
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FIG. 2. The ratio of the total width 8 to t and the ra- 
diative width Srad of TM;,i WG resonances is shown as a 
function of the imaginary part ni of the refractive index for 
two values of the real part ur of the refractive index and 
various values of /. The values 1 = 70, 60, 50 correspond to 
i?/A~8.4, 7.3, 6.2, respectively, and the values 1 = 22, 20, 15 
correspond to J?/A~1.8, 1.6, 1.3, respectively. 

Introducing the quality factor Q = £1/(2(5), the ab- 
sorption-assisted part Qabs — ^/(2<5 a bs) is derived from 
Eqs. 0-© to be 

Qabs -2^ + 2^(nl-ir/^ + ° {V > (20) 



In order to obtain the exact quality factor in de- 
pendence on the resonance frequency, the permittiv- 
ity as a function of frequency must be known and the 
calculations should be based on the original equations 
( |j~0| ) and (|ll]). In particular, the quality factor Q ra d, 
which accounts for the line broadening associated with 
the input-output coupling, can be obtained by setting 
the imaginary part of the permittivity equal to zero, 
so that Q — Q r ad- Recalling Eq. (|l^), we may write 
1/Q = l/Qtot = 1/Qrad + 1/Qabs, from which the quality 
factor Qabs, which accounts for the line broadening due 
to material absorption, can then be calculated. Typical 
results obtained for TM WG waves on the basis of the 
model permittivity in Eq. (Q) are shown in Fig. 

From Fig. ||(a) it is seen that (for chosen radius of the 
microsphere) Q a bs decreases with increasing frequency 
(i.e., increasing nj ~ e/), whereas Q ra d increases with 
increasing frequency (i.e., increasing ur ~ \/£r)- Suffi- 
ciently below the band gap where tir is small and the 
radiative losses dominate, Qtot follows Q ra d- With in- 
creasing frequency the absorption losses increase and can 
eventually dominate the radiative losses, so that now Qtot 
follows Q abs . 

In Fig. ||(b), the values of Qabs calculated from the 
exact equation ( ji~l| ) are compared with the approximate 
values according to Eq. ( |20| ) . It is seen that the neglect 
of dispersion in Eq. ( p0[ ) leads to some underestimation 
of the absorption-assisted quality factor. The difference 
between the exact and the approximate values increase 
with frequency, because of the increasing dispersion. We 
have also calculated Qabs from the approximate equation 
( lilt ) [together with Eq. (|i|)] and found a good agreement 
with the exact result. 

Figure ^(c) illustrates the influence on the quality fac- 
tor of the radius of microsphere and the strength of ab- 
sorption. As expected, for lower frequencies where Qtot 
follows Q ra d, the quality is improved if the radius of the 
sphere is increased. For higher frequencies where Qtot 
follows Qabs , the quality becomes independent of the ra- 
dius. Note that with increasing radius the (frequency) 
distance between neighboring resonances decreases. 



[P = n R (l/n R ) for TE (TM) waves]. When the absorp- 
tion losses dominate the radiative ones, it is frequently 
assumed [^|Jl9|l that the quality factor is simply given by 
the ratio of the (plane-wave) absorption length and the 
wavelength, i.e., Q = nn/(2ni), which is just the first 
term on the right-hand side in Eq. (|o|). As long as the 
values of v and/or ur are large, which is the case for WG 
waves, the second term on the right-hand side in Eq. j20| ) 
is typically less than a few percent, and the first term is 
indeed the leading term. Nevertheless, Eq. ( p0| ) may be- 
come significantly wrong, because dispersion is fully dis- 
regarded. Moreover, the radiative losses can dominate 
material absorption. 



B. Resonances inside the band gap 

Let us now look for (transverse) field resonances inside 
the band gap, which is a strictly forbidden zone for non- 
absorbing bulk material. Inside the gap, we may assume 
that the real part of the refractive index is much smaller 
than the imaginary part, hr^uj) <^ni(uj). Using Bessel- 
function expansion (Appendix B), it can be shown that 
there are no TE resonances [Eq. ( |i"o| ) has no solutions], 
and the complex roots of Eq. (^) , which determine the 
TM resonances, can formally be given by 
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FIG. 3. Quality factors of TMz : i WG field resonances in a 
dielectric microsphere of radius R and complex permittivity 
e(w) according to Eq. (Q) (lop = 0.5 lot)- (a) Exact values of 
Qrad, Qabs, and Qtot for fl=10A T , 7 = 10 _6 a;T, 48<Z<78. 
(b) Exact values of Qabs (+) and approximate values (A) ac- 
cording to Eq. ( pcf ) for the same parameters as in (a), (c) 
Exact values of Qtot for 7 /lot = 10~ 6 (o) and j/lot = 10~ 4 
(A); (1) ,R=10A T , 48<Z<78; (2) R=2X T , 10 < Z < 22. 

e(w ; ) 2 + e(^;) + l 
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I/W1+ 
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e(wj) 2e(w z ) A /-e(w i )-l 



OK 1 ), 

(21) 



where the condition 
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FIG. 4. Quality factors of TM; SG field resonances in a 
dielectric microsphere of radius 7? and complex permittivity 
e(io) according to Eq. (|l|) (up — 0.5 lot)- (a) Exact values of 
Qrad, Qabs, and Q to t for i? = 10A T , j = W a u T , 70<Z<120. 
(b) Exact values of Q a b s (+) and approximate values (A) ac- 
cording to Eq. ( |24| ) for the same parameters as in (a), (c) 
Exact values of Qtot for j/ut — W -6 (°) and j/ut = 10 -4 
(A); (1) R = WX T , 70</<86; (2) R = 2 X T 15 < I < 55. The 
gap is from o;t to ~ 1.12 ut [cf. Fig.hJ . According to 
Eq. d23|), resonances exist for f2; < 1.0607 ojt. 



£i?(w;) < -1 



(22) 



has been required to be satisfied. Note that for the 
Drude-Lorentz model (0), condition ( p2| ) leads to 
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(23) 



The total linewidth can be determined according to 
Eq. (17), which is still valid here. If we ignore the depen- 
dence on frequency of the permittivity, then the expres- 
sion on the right-hand side of Eq. ( pl| ) can be regarded 
as a function of the refractive index. The contribution 
to the linewidth of the material absorption can then be 
calculated by first-order Taylor expansion at ni of this 
function in close analogy to Eq. ( |l6| ) (with the roles of tir 
and ni being exchanged) . In terms of the quality factor 
Qabs = fi/ (2<5 a f«) the result reads 



Qa 



bs 



nj(n^-l) 
2n R 



+ 0(0- 



(24) 



Since now tir is proportional to 7, Q a b s is again propor- 
tional to 7 _1 . 

It is well known that when the condition (|2|) is sat- 
isfied, then surface-guided (SG) waves can be excited - 
waves that are bound to the interface and whose am- 
plitudes are damped into either of the neighboring me- 
dia [ 2^] . Typical examples are surface phonon polaritons 
for dielectrics and surface plasmon polaritons for metals. 
Obviously, the TM resonances as determined by Eq. ( |2l| ) 
correspond to SG waves of a sphere. Note that for the 
SG waves, in contrast to the WG waves, each angular 
momentum number I is associated with only one wave. 

Examples of the frequencies and quality factors of SG 
waves are plotted in Fig. || From a comparison with 
Fig. |^ it is seen that the quality factors of the SG waves 
are comparable with those of the WG waves and also 
behave in a similar way as the latter. So, the radiative 
losses are the dominant ones for lower-order resonances, 
whereas for higher-order resonances the losses essentially 
arise from material absorption [Figs. |](a) and (c)], and 
the radiative losses can be reduced by increasing the ra- 
dius of the sphere [Fig. ||(c)]. Further, a neglect of dis- 
persion again leads to some overestimation of material 
absorption [Fig. §(b)]. 



C. Ground-state fluctuations 



When an excited atom is in free space, then its cou- 
pling to the vacuum-field fluctuation gives rise to spon- 
taneous decay. In the presence of dielectric bodies, the 
fluctuation of the electromagnetic field with respect to 
the ground-state of the combined system that consists of 
the radiation field and the dielectric matter must be con- 
sidered. The electric-field correlation in the ground state 
can be characterized by the correlation function 
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FIG. 5. Spatial variation of the ground-state fluctuations 
P rr {r,uj) (in arbitrary units) inside a dielectric microsphere of 
radius R = 2Xt and complex permittivity e(u) according to 
Eq. (0) [lu p = 0.5ljt; = 10" 4 (solid line), 10" 6 (dashed 

line), and 1(T 8 (dotted line)], (a) w = 0.94042 w T , TM 16 ,i WG 
field resonance, (b) u> — 1.02811 lot, TM16 SG field resonance, 
(c) w = 1.05339 u T , TM 30 SG field resonance. 
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FIG. 6. The same as in Fig. |f| but for outside the sphere. 

where E^ir, ui) is the electric- field operator in the fre- 
quency domain and Gy(r,r',u) is the classical Green 
tensor of the dielectric-matter formation of given complex 
permittivity e(r, ui) that satisfies the Kramers-Kronig re- 
lations. For equal spatial arguments, the power spectrum 
of the ground-state fluctuation can be obtained according 
to the relation 



(0|^(r,w)4(iV)|0) = P y (r,«)J(w - w') 
which together with Eq. ( |25| ) reveals that 

hu 2 



Im Gij(r,r,Lu). 



(26) 



(27) 



From the linearity of the Maxwell equations it follows 
that the Green tenso r fo r a dielectric body can be de- 
composed as [cf. Eq. (Al)] 



Ga (r, r', «) = G°, (r, r', u) + G« (r, r', w), 



(28) 



where G°- is the Green tensor for either the vacuum (out- 
side the body) or the bulk material (inside the body) , and 
the scattering term Gfj ensures the correct boundary con- 
ditions (at the surface of discontinuity). Inside the body 
ImG'j becomes singular for r— >r', and some rcgular- 
ization (e.g., averaging over a small volume) is required. 
Since we are only interested in the spatial variation of the 
ground-state fluctuation, we may disregard the (space- 
independent) contribution that arises from Im G? . 

The spatial variation of the radial diagonal component 
P rr inside and outside the microsphere under consider- 
ation is illustrated in Figs. |5| and respectively. Note 
that P rr only refers to the TM field noise [cf. Eqs. (A3) 
- (A7)]- It is seen that the fluctuation associated with 



the WG-type field [Figs. ||(a) and ||(a)] essentially con- 
centrates inside the sphere, whereas the fluctuation as- 
sociated with the SG-type field [Figs. ||(b,c) and ||(b,c)] 
concentrates in the vicinity of the surface of the sphere. 
In both cases, a singularity is observed at the surface. In 
Figs. ||(a) and ||(a), we have restricted our attention to a 
single-maximum WG field, i.e., i=l. Otherwise i (> 1) 
maxima would be observed inside the sphere. 

The strong enhancement of the fluctuation which leads 
to the singularity at the surface comes from absorption 
and gives rise to nonradiative energ y transfer from the 



[II B). Clearly, for suffi- 



atom to the medium (see Sec. 
ciently small absorption the divergence becomes mean- 
ingless, because the required (small) spatial resolution 
would be illicit. From Figs. ^|(a) and ^(a) it is seen that 
material absorption reduces the fluctuation of the WG 
field inside and outside the sphere, and Figs. |5|(b,c) and 
^(b,c) show that absorption changes the area on either 
side of the surface over which the SG field fluctuation 
extends. As expected, the effects are less pronounced 
for low-order resonances where the radiative losses are 
the dominant ones [compare Figs. |^(b) and (c) with 
Figs. ||(b) and (c) respectively]. 



III. SPONTANEOUS DECAY 



A. Basic formulae 



In the weak-coupling regime, an excited atom near a 
a dispersing and absorbing body decays exponentially, 
where the decay rate is given by 



tier, 



(29) 



7 



(k A — w A /c; [a, atomic transition dipole moment), and 
the contribution of the body to the shift of the transition 
frequency reads 



irnco 



duj 



; 2 lmGfj{r A ,r A ,L 
: 2 uj — U) A 



(30) 



which can be rewritten, on using the Kramers-Kronig re- 
lations, as 



Sua = 



he 



V 



du> 



ReGfj{r A ,Y A ,ujA) 
to 2 lmGfAr A ,r A ,uj)' 



■<jj A 



(31) 



(for details, see |I^,|T^]). It is not difficult to see that 
in Eq. ( pl| ) the second term, which is not sensitive to 
the atomic transition frequency, is small compared to the 
first one and can therefore be neglected. Note that the 
vacuum Lamb shift (see, e.g., Mj) may be thought of as 
being already included in the atomic transition frequency. 

The intensity of the spontaneously emitted light regis- 
tered by a p ointlike photodetector at position r and time 
t reads 



k A fij 



ire 



dt' 



C u (t') 



did Im Gij (r , , uj) e 



-i{uj—ujA){t—t') 



(32) 



where C u (t) is the probability amplitude of finding the 
atom in the upper state. Note that Eq. d32| ) is valid for 
an arbitrary coupling regime. In particular, in the weak 
coupling regime, where the Markov approximation ap- 
plies, C u (t') can be taken at t' = t and put in front of the 
time integral in Eq. (|32|), with C u (t) being simply the 
exponential 



C u {t) = e (-A/2+iSu A )t 

Equation ( |32|) thus simplifies to 

I{v,t)^\F{v,v Al uj A )\ 2 e- A \ 

where 



(33) 



(34) 



Fi{v, r A ,w A ) 



V 



du 



Im G i:j (y,T A ,uj) 

LO + UJ A 



Gij(T,T A ,CJ A ) 



(35) 



Since the second term on the right-hand side of Eq. ( |35| ) 
is small compared to the first one, it can be omitted, 
and the spatial distribution of the emitted light (emis- 
sion pattern) can be given by, on disregarding transit 
time delay, 



\F(r,r A ,uj A ) 



E 



^ A ^3 



G y -(r, r A ,uj A ) 



(36) 



Material absorption gives rise to nonradiative decay. 
Obviously, the rate (|2|) as the total decay rate describes 
both radiative decay and nonradiative decay. The frac- 
tion of emitted radiation energy can be obtained by in- 
tegration of I(r, t) with respect to time and integration 
over the surface of a sphere whose radius is much larger 
than the extension of the system, 



W = 2ce 



dt 



d<f) / dOp 2 sinO I(r,t) (37) 



(p = r — r^). The ratio W/W , where W = Tiuj a is the 
emitted energy in free space, then gives us a measure of 
the emitted energy, and accordingly, 1 — W/Wq measures 
the energy absorbed by the body. 



B. Spontaneous decay rate 



Usmg Eqs. (|AfJ) - (|A3J) and QA§) - (|A7j), the decay 
rate (|29| ) for a (with respect to the microsphere) radially 
oriented transition dipole moment can be given by 

f °° 
A x = A h + §£)l(J + l)(2J+l) 



x Re 



(38) 



and for a tangential dipole it reads 

{oo 
l + f^(2I + l) 
i=i 

Bf\oJ A ) (hf\k A r A ) 

[k A r A h\ 1 \k A r A ''~ 1 ' ■ 2 



x Re 



k A r A 



(39) 



where the prime indicates the derivative with respect to 
k A r A , and 



A = 



(40) 



is the rate of spontaneous emission in free space. Note 
that a radially oriented transition dipole moment only 
couples to TM waves, whereas a tangentially oriented 
dipole moment couples to both TM and TE waves. Equa- 
tions ([38]) and (|39| ) generalize the results obtained for 
nondispersing and nonabsorbing matter whose resonance 
frequencies are far from the atomic transition frequencies, 
i.e., e = eji > 1, [ |lO|]TT| to arbitrary Kramers-Kronig con- 
sistent matter, without placing restrictions to the transi- 
tion frequency. 
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FIG. 7. Decay rate versus atomic transition frequency 
for a radially oscillating dipole near a microsphere of ra- 
dius R = 2\t and complex permittivity e(ui) according to 
Eq. ([j} (ljp = 0.5lut, -y/uj T = W 4 )- (a) Ar = 0.02A T ; in- 
set: 7/aj T =lCr 4 (solid line), 10" 5 (dashed line) , 10" 6 (dotted 
line), (b) Ar = 0.1 At. 



atom from the surface of the microsphere must not be 
smaller than interatomic distances (otherwise the con- 
cept of macroscopic electrodynamics fails), the diver- 
gence at the surface (Ar = 0) is not observed. 

The dependence of the decay rate (38) on the radiat- 
ion-field structure, which can be observed for not too 
small (large) values of Ar (ej), is illustrated in Fig. |t] for 
a radially oriented transition dipole moment. Since the 
decay rate is proportional to the spectral density of final 
states, its dependence on the transition frequency mimics 
the excitation spectrum of the sphere-assisted radiation 
field. As a result, both the WG and SG field resonances 
can strongly enhance the spontaneous decay. The en- 
hancement decreases with increasing distance between 
the atom and the sphere [compare Figs. ||(a) an(b)]. 
When material absorption increases, then the resonance 
lines are broadened at the expense of the heights and 
the enhancement is accordingly reduced [see the inset in 
Fig. @(a)] 

Figure M reveals that the SG field resonances, which 
are the more intense ones in general, can give rise to 
a much stronger enhancement of the spontaneous decay 
than the WG field resonances. In particular, with in- 
creasing angular-momentum number the lines of the SG 
field resonances strongly overlap and huge enhancement 
factors can be observed for transition frequencies inside 
the band gap [e.g., factors of the order of magnitude 
of 10 4 for the parameters chosen in Fig. 0(a)]. When 
the distance between the atom and the sphere increases, 
then the atom rapidly decouples from that part of the 
field. Thus, the huge enhancement of spontaneous decay 
rapidly reduces and the interval in which inhibition of 
spontaneous decay is typically observed, extends accord- 
ingly [see Fig. @(b)]. 



When the atom is very close to the microsphere, 
Eqs. ( p8| ) and ( |39| ) simplify to (see Appendix |C]) 



3A ei(u A ) 



1 



4fci \e(u A ) 



1| 2 (Ar) s 



■0(1) 



(41) 



and 



A \\ = 3A c^a) 



8k 3 A |e(u; A ) 



1 



1 



(Ar) 



2R 2 Ar 



O(l), (42) 



with Ar = ta~R being the distance between the atom 
and the surface of the microsphere. The terms ~ (Ar) -3 
and oj (Ar) -1 result from the TM near- field coupling. 
Since they are proportional to e/, they describe non- 
radiative decay, i.e., energy transfer from the atom 
to the medium, and reflect the strong enhancement 
of t he e lectromagnetic-field fluctuation as Ar — > (see 
Sec. PQ). In particular, the terms ~ (Ar)~ 3 in Eqs. ( J4l] ) 
and (W2) are exactly the same as in the case when the 
atom is close to a planar interface |22|,^3|]. Obviously, 
nonradiative decay does not respond sensitively to the 
actual radiation-field structure. Since the distance of the 



C. Lamb shift 



Substituting expressions flA3| ) and QA5| ) - (A7) into 
Eq. (|3|), we find that the shift of the atomic transition 
frequency, which is caused by the presence of the micro- 
sphere, is 



3A 



£i(i + l)(2*+l) 



Im B?{w A ) 



hP(k A r A ) 



-1 2 



k A T A 



(43) 



for a radially oriented transition dipole moment, and 



*4 = -^£(2i + i) 



i=i 



x Imi Bf{tUA 



[hf\k A r A ) 
[kATAh^jkArA)}' 

kATA 



(44) 
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for a tangentially oscillating dipole. In particular, when 
the distance between the atom and the sphere is very 
small, then Eqs. ( fl3| ) and ( f44| ) simplify to [in close anal- 
ogy to Eqs. © and @] 



Suji 



3A a \e(lo a )\ 2 - 1 1 
16fc3 |e( Wj4 ) + l| 2 (Ar) s 



O(l) 



and 

sJ 



3A \e(u A )\ 2 -l 
32k A \e{u A ) + l\ 2 



(Ar) 3 2i? 2 Ar 



(45) 

0(1). 

(46) 



The terms ~ (Ar)~ 3 and ~ (Ar) -1 again result from the 
TM near-field coupling, and the leading terms agree with 
those obtained for a planar interface p3[. Note that in 
contrast to the decay rate, the Lamb shift diverges for 
Ar — > even when ej = 0. 

The dependence of the frequency shift ([l3|) on the rad- 
iation-field structure for not too small values of Ar is 
illustrated in Fig. ||. It is seen that each field resonance 
can give rise to a noticeable frequency shift in the very 
vicinity of the resonance frequency; transition frequen- 
cies that are lower (higher) than the resonance frequency 
are shifted to lower (higher) frequencies. In close anal- 
ogy to the behavior of the decay rate, the shift is more 
pronounced for SG field resonances than for WG field 
resonances and can be huge for large angular momen- 
tum numbers when the lines of the SG field resonances 
strongly overlap. 

The behavior of the Lamb shift as shown in Fi g. js| (b) 
can already be seen in the single-resonance limit pjF , If 
the atomic transition frequency u>a is close to a reso- 
nance frequency Q of the microsphere, contribution from 
other resonances may be ignored in a first approximation. 
Regarding the resonance line as being a Lorentzian and 
using contour integration, we obtain from Eq. (BQ) 



5u>a 



A{n)8 2 

47T 

A(Q)6 



1 



1 



uj — uja (lu — Q,) 2 + 5 2 



A 



S 2 



(47) 



where A(f2) is the decay rate as given by Eq. ( p9|) (with 
in place of to a), and A — uja — CI. In particular, Eq. ( f47j ) 
indicates that the frequency shift peaks at half maximum 
on both sides of the resonance line. With increasing ma- 
terial absorption, the linewidth S increases while A(Cl) 
decreases and thus the absolute values of the frequency 
shift are reduced, the distance between the maximum 
and the minimum being somewhat increased. With de- 
creasing distance between the atom and the microsphere 
near-field effects become important and Eq. ( f47j ) fails, as 
it can be seen from a comparison of Figs. |](a) and (b). 
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FIG. 8. Lamb shift versus atomic transition frequency 
for a radially oscillating dipole near a microsphere of ra- 
dius R = 2\t and complex permittivity e(u)) according to 
Eq. (|l|) (up =0.5 lu t , 7/w T = 10~ 4 ). (a) Ar = 0.02A T . 
(b) Ar = 0.1A T ; inset: j/lj t =1(T 4 (solid line), 1(T 5 (dashed 
line), 10" 6 (dotted line). 



D. Emitted-light intensity 

1. Spatial distribution 



Using Eqs. ( |aT| ) - ( |A7|) , the function F(r,rj,w A ), 
Eq. (|35|), which determines, according to Eq. (|34|), the 
spatial distribution of the emitted light, can be given by 
io A = <t>A = 0,r A <r) 



kATA 



1 = 1 



ji(k A rA) + Bj\iUA)h' i 1 \kArA) 



k A r 



[k A rh^\k A r) 
kir 



sin 6»P/(cos( 



(48) 



for a radially oriented transition dipole moment, and 
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-11/ \ k 3 Af i^(2l + l) 

F " (r ' r -^^g7(rrTy 

x| e r cos 4>B? 1(1+1) fe ' ] } kA ^ sin0P/(cos0) 
I Mr 



- eg COS ( 



B ; M /i ; (1) (fc A r)P/(cos0) 



B - r [^fM]'p [(cos , 



e„i sin < 



S ; M /i ; (1) (fc j4 r)P i (cos( 



5Ar [^rftf'tM]' D , 



k A r 



P/(cos( 



(49) 



for a tangentially oriented dipolc in the xz-plane. Here, 
the abbreviating notations 



B? = 



— \[kArA3i{kArA)]' 
\r a y 



^ATA 



+ B?(u A )[k A r A hf ) (k A r A )]'} 1 



= j l {kArA) + B^A)h\ 1 \k A rA), 



(50) 



(51) 



P(cos0) = 1(1 + l)P/(cos0) - cos 0P{ (cos 0) (52) 

have been introduced. Examples of the far-field emission 
pattern of a radially oriented transition dipole moment, 
iF^r, r^,o;^)| 2 , and a tangentially oriented transition 
dipole moment, |FH(r,r^,w^)| 2 , are plotted in Figs. ^ 
and [lC] respectively. 

Let us first restrict our attention to a radially oriented 
transition dipole moment. In this case, the far field is es- 
sentially determined by Fg- , as an inspection of Eq. ( |4S| ) 
reveals. When the atomic transition frequency coincides 
with the frequency of a WG wave of angular momen- 
tum number I far from the band gap [Fig. 0(a)], then the 
corresponding Z-term in the series (Eq) obviously yields 
the leading contribution to the emitted radiation, whose 
angular distribution is significantly determined by the 
term ~sin0P/(cos0). Since Pi(cos9) is a polynomial 
with I real, single roots in the interval < 9 < n [p5| , 
sin 8 P/(cos 9) must have I extrema within this interval. 
Thus, the emission pattern has I lobes in, say, the yz- 
plane, i.e., I cone-shaped peaks around the z-axis, be- 
cause of symmetry reasons. The lobes near 9 = and 
9 = 7r are the most dominant ones in general, because of 

Pel 
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sin0P/(cos0) ~ (sin0)- 1/2 + 0(1^) 



(53) 



(0 < 9 < 7r). Note that the superposition of the leading 
term with the remaining terms in the series (|48| ) gives rise 
to some asymmetry with respect to the plane 9 = it/ 2. 



FIG. 9. Polar diagrams of the normalized far-field emis- 
sion pattern |F x (r, ta, oja)\ 2 /(k A ^/4:TYeo) 2 of a radially os- 
cillating dipole near a microsphere of radius R=2Xt and 
complex permittivity e(u) according to Eq. (|l]) (top = 0.5 ujt, 
^/ljt = W~ 4 , Ar = 0.02X T , r = 20X T ). uj A /oJT = {a) 0.94042, 
(b) 0.999, (c) 1.02811, and (d) 1.06. 
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FIG. 10. The same as in Fig. but for a tangentially 
oscillating dipole. 

When the atomic transition frequency approaches the 
band gap (but is still outside it), a strikingly different 
behavior is observed [Fig. ^(b)]. The emission pattern 
changes to a two-lobe structure similar to that observed 
in free space, but bent away from the microsphere sur- 
face, the emission intensity being very small. Since near 



the band gap absorption losses dominate, a photon that 
is resonantly emitted is almost certainly absorbed and 
does not contribute to the far field in general. If the pho- 
ton is emitted in a lower-order WG wave where radiative 
losses dominate, it has a bigger chance to escape. The 
superposition of all these weak (off-resonant) contribu- 
tions just form the two-lobe emission pattern observed, 
as it can also be seen by careful inspection of the series 

(ED- 

When the atomic transition frequency is inside the 
band gap and coincides with the frequency of a SG wave 
of low order such that the radiative losses dominate, then 
the emission pattern resembles that observed for reso- 
nant interaction with a low-order WG wave [compare 
Figs. ||(a) and (c)]. With increasing transition frequency 
the absorption losses become substantial and eventually 
change the emission pattern in a quite similar way as do 
below the band gap [compare Figs. ||(b) and (d)]. Obvi- 
ously, the respective explanations are similar in the two 
cases. 

For a tangentially oriented transition dipole moment 
the situation is not lucid in general. Let us therefore re- 
strict our attention to the far field in the xz-plane, i.e., 
cf> — in Eq. (|49|). In this case, the main contribution 

— II 

to the far field comes from Fg. The interpretation of 

the plots in Fig. |l^ is quite similar to that of the plots 
in Fig. ^. In particular, when the atomic transition fre- 
quency coincides with a WG field resonance that mainly 
suffers from radiative losses, then the I term in the series 
( Jl9| ) that corresponds to the order of the WG wave is the 
leading one, and the main contribution to it stems from 
the term ~ 1(1 + l)P/(cos 0). It gives rise to I — 1 lobes in 
the interval < 8 < tt, and two lobes at 9 = and 8 = it, 
which are the most pronounced ones [Fig. |lO](a)]. These 
maxima are approximately located at the positions of the 
minima of the far field of the radially oscillating dipole. 
Hence, if the transition dipole moment has a radial and 
a tangential component, a smoothed superimposed field 
is observed. 



2. Radiative versus nonradiative decay 

Since both the imaginary part of vacuum Green ten- 
sor G?„ and the scattering term are transverse, the 

I— H l J 

decay rate Q29| ) results from the coupling of the atom to 
the transverse part of the electromagnetic field. Never- 
theless, the decay of the excited atomic state must not 
necessarily be accompanied by the emission of a real pho- 
ton, but instead a matter quantum can be created, be- 
cause of material absorption. To compare the two decay 
channels, we have calculated, according to Eq. (|37|), the 
fraction W/Wq of the atomic (transition) energy that is 



c (tra: 

irradiated. Using Eqs. ( |34]), (|37| ), a nd Q4S| ), we derive, 
on recalling the relations (A13) and (A14), for a radially 



oriented transition dipole moment 
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W 3A, 



W 2A A 
1 



£l(Z + l)(2I + l) 



1=1 



(k A r A ) 2 



\ji{k A r A )+Bf(u JA )h ( l i \k A r A ) . (54) 



Recall that W/Wq = 1 implies fully radiative decay, while 
W/Wq = fully nonradiative one. 

The dependence of the ratio W/Wq on the atomic tran- 
sition frequency is illustrated in Fig. [ll]. The minima at 
the WG field resonance frequencies indicate that the non- 
radiative decay is enhanced relative to the radiative one. 
Obviously, photons at these frequencies are captured in- 
side the microsphere for some time, and hence the prob- 
ability of photon absorption is increased. For transition 
frequencies inside the band gap, two regions can be dis- 
tinguished. In the low-frequency region where low-order 
SG waves are typically excited radiative decay dominates. 
Here, the light penetration depth into the sphere is small 
and the probability of a photon being absorbed is small 
as well. With increasing atomic transition frequency the 
penetration depth increases and the chance of a photon 



to escape drastically diminishes. As a result, nonradia- 
tive decay dominates. Clearly, the strength of the effect 
decreases with decreasing material absorption [compare 
Figs. [y](a) and (c) with Figs. |TT|(b) and (d) respectively]. 

From Figs. pd|(a) and (b) two well pronounced min- 
ima of the totally emitted light energy, i.e., noticeable 
maxima of the energy transfer to the matter, are seen for 
transition frequencies inside the band gap. The first min- 
imum results from the overlapping high-order SG waves 
that mainly underly absorption losses. The second one 
is observed at the longitudinal resonance frequency of 
the medium. It can be attributed to the atomic near- 
field interaction with the longitudinal component of the 
medium-assisted electromagnetic field, the strength of 
the longitudinal field resonance being proportional to e/. 
Hence, the dip at the longitudinal frequency of the emit- 
ted radiation energy reduces with decreasing material ab- 
sorption [compare Fig. |ll](a) with Fig. p!l|(b)], and it dis- 
appears when the atom is moved sufficiently away from 
the surface [compare Figs. |ll](a) and (b) with Figs. [fl](c) 
and (d) respectively]. 



13 



0.01 



a 

CD 



2 

"3 



0.005 







1 





0.05 



A t 



0.1 



FIG. 12. Exact [Eq. (p2|), solid lines] and approxi- 
mate [Eq. (^), dashed lines] time evolution of the far-field 
intensity I(r, t)/(k A fi/4ireo) 2 at a fixed point of observa- 
tion for a radially oscillating transition dipole moment 
(A /lo t = W~ 7 , R = 2\ T , Ar = 0.02A T , r = 20A T , = 3, 



lop — 0.5ujt; 7=10" 
wave with Q~ 10 3 ). 
with Q~10 4 ). 



4 w T ). (a) lua = 0.91779 lot (TM u ,i WG 
(b) lj a = 0.94042 cot (TMie.i WG wave 



3. Temporal evolution 

Throughout this paper we have restricted our attention 
to the weak-coupling regime where the excited atomic 
state decays exponentially, Eq. (|33|). When retardation 
is disregarded, then the intensity of the emitted light (at 
some chosen space point) simply decreases exponentially, 
Eq. (34). To study the effect of retardation, we have also 
performed the frequency integral in the exact equation 
( |32| ) numerically. 

Typical examples of the temporal evolution of the far- 
field intensity are shown in Fig. 12 for a radially oriented 
transition dipole moment in the case when the atomic 
transition frequency coincides with the frequency of a 
WG wave. Whereas the long-time behavior of the inten- 
sity of the emitted light is, with little error, exponential, 
the short-time behavior (on a time scale given by the 
atomic decay time) sensitively depends on the quality fac- 
tor. The observed delay between the upper-state atomic 
population and the intensity of the emitted light can be 
quite large for a high-Q microsphere, because the time 
that a photon spends in the sphere increases with the 
Q value. Further, in the short-time domain some kink- 
like fine structure is observed, which obviously reflects 
the different arrival times associated with multiple re- 
flections. 

The results in Fig. |lj refer to a fixed space point. 
Figure [l3] illustrates the transient behavior of the spa- 
tial distribution of the emitted light. In particular, it is 
seen that some time is necessary to build up the sphere- 
assisted spatial distribution of the emitted light which is 
typically observed for longer times when the approxima- 
tion (|34|) applies. 
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FIG. 13. Exact [Eq. (g^), solid lines] and approximate 
[Eq. (|34|), dashed lines] angular distribution of the far-field 
intensity 7(r, t) / {k\[i/ Aitea) 2 at different times for a radially 
oscillating transition dipole moment [for the parameters, see 
Fig. 0, curves (b)]. 



IV. METALLIC MICROSPHERE 

Setting in Eq. @) lot = 0, we obtain (within the Drude- 
Lorentz model) the permittivity of a metal. Hence, the 
results derived for the band gap of a dielectric micro- 
sphere also applies, for appropriately chosen values of 
LOp and 7, to a metallic sphere. The examples plotted in 
Figs. [14] and [15] refer to silver (u P = 1.32 x 10 16 s _1 and 
1/7=1.45 x 10- 14 s §f|). 

Figure ^ shows the positions and quality factors of SG 
field resonances for different microsphere radii. The be- 
havior is quite similar to that observed for a dielectric mi- 
crosphere [cf. Fig. Q]. The radiative losses decrease with 
increasing radius of the sphere, while the losses due to 
material absorption are less sensitive to the microsphere 
size [compare Figs. |l4](a) and (b)]. It is further seen that 
Qrad increases with the angular momentum number of 
the SG wave, while Q a bs slightly decreases. 
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FIG. 14. Quality factors Q rad (A), Q abs (+), and 
Qtot (°) of TM; SG field resonances in a metallic mi- 
crosphere of complex e(cu) [Eq. ([j]), lut~0, j/ujp — 0.005]. 
(a) R=10\ P , 30</<60. (b)i?=5A P , 16</<60. All res- 
onances obey, in accordance with the condition (|2^), the re- 
lation uo/ujp < l/\/2^0.71. 
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FIG. 15. (a) Decay rate of an atom near a a metallic 
microsphere of complex permittivity e(u) [Eq. (^), uit~0, 
y/ujp — 0.005] as a function of the transition frequency 
for a radially oriented transition dipole moment (R — 5Xp, 
Ar — 0.1 Xp). (b) Polar diagram of the normalized far-field 
emission pattern |F x (r, r^, lua)\ 2 /(k A (i/iireo) 2 for r = 50Ap 
and uja/ijJp =0.5026, the other parameters being the same as 
in (a). 



As a result, the radiative losses again dominate for 
low-order resonances and the absorption losses for high- 
order ones. Since absorption tends to be larger in metals 
than in dielectrics, the dominance of material absorption 
can already set in at lower-order resonances. Note that 
even for metals the relationship (||) typically still holds, 
so that Eqs. @ - (§|) apply. 

The dependence of the decay rate on the transition fre- 
quency of an excited atom placed near a metallic micro- 
sphere is illustrated in Fig. |l5](a) for a radially oriented 
transition dipole moment. An example of the emission 
pattern for the case when the atomic transition frequency 
coincides with the frequency of a SG wave is shown in 
Fig. |l5|(b) . When the radius of the microsphere becomes 
too small, then SG waves cannot be excited. In par- 
ticular, in ]2^] it was assumed that R <^ Ap, and thus 
the resonances shown in Figs. |IJ and [l5](a) could not be 
found. It is worth noting that, in contrast to dielectric 
matter, a large absorption in metals can substantially 



enhance the near-surface divergence of the decay rate, 
[Eqs. ( fll| ) and (fl2|)], which is in agreement with experi- 
mental observations of the fluorescence from a thin layer 
of optically excited organic-dye molecules that were sep- 
arated from a planar metal surface by a dielectric layer 
of known thickness p{| . 



V. CONCLUSIONS 

We have applied the recently developed formalism |l3| ] 
to the problem of spontaneous decay of an excited atom 
near a dispersing and absorbing microsphere. Basing the 
calculations on a complex permittivity of Drude-Lorentz 
type, which satisfies the Kramers-Kronig relations, we 
have been able to study the dependence of the decay 
rate on the transition frequency for arbitrary frequen- 
cies. We have shown that the decay can be substantially 



15 



enhanced when the transition frequency is tuned to either 
a WG field resonance below the band gap (for a dielectric 
sphere) or a SG field resonance inside the band gap (for 
a dielectric sphere or a metallic sphere). 

Whereas for both low-order WG field resonances and 
low-order SG field resonances radiative losses dominate, 
high-order field resonances mainly suffer from material 
absorption. Accordingly, spontaneous decay changes 
from being mainly radiative to being mainly nonradia- 
tive, when the transition frequency (tuned to a field reso- 
nance) in the respective frequency interval increases. We 
have further shown that in the presence of strong ma- 
terial absorption the decay rate drastically raises as the 
atom approaches the surface of the microsphere, because 
of near-field assisted energy transfer from the atom to 
the medium. Thus, the effect is typically observed for 
metals. 

When radiative losses dominate, the emission pattern 
is highly structured, and a substantial fraction of the 
light is emitted backward and forward within small po- 
lar angles with respect to the tie line between the atom 
and the center of the sphere. With increasing absorption 
this directional characteristic is lost. When absorption 
losses dominate, the (weakened) emission pattern takes 
a form that is typical of reflection at a mirror. 

In the paper we have restricted our attention to the 
weak-coupling regime, assuming that the excited atomic 
state decays exponentially. Obviously, when the atomic 
transition frequency coincides with a resonance frequency 
of the cavity-assisted field, the strong-coupling regime 
may be realized. In particular, SG waves seem to be best 
suited for that regime, because of the noticeable enhance- 
ment of spontaneous emission. The calculations could be 
performed in a similar way as in Ref. |l3fl. 
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APPENDIX A: THE GREEN TENSOR 

The Green tensor of a microsphere of radius R (reg- 
ion 2) embedded in vacuum (region 1) can be decom- 
posed into two parts p0 |, 

G(r, r', u) = G< s >(r, r', u)8 fa + G^>(r, r', w), (Al) 

where G^ (r, r', co) represents the contribution of the di- 
rect waves from the source in an unbounded space, and 
GU s )(r, r',w) is the scattering part that describes the 
contribution of the multiple reflection (/ = s) and trans- 
mission (/ ^ s) waves (/ and s, respectively, refer to the 



regions where are the field and source points r and rO. 
In particular, G^ , G^ n \ and G^ can be given by @] 



G (s) (r,r', 



d(r-r') 



(2 — Som) 



ik s \ ~\ 21 + 1 (I — to)! 

: [Mf lm (r,k s )M llm (r',k s ) 
+ ^P lm (r,k s )N llm (r',k s ) 



if r>r', and G (s \r, r', w) = G^( r ': r,u>) if r<r', 
G( n )(r,r» 

4 V^V^ V" 2/ + 1 (I — m)\ 



(A2) 



1 = 1 771=0 



1(1 + 1) (l + m)\ 



(2— So m ) 



[Br^)Mf lm (r,ki)Mf lm (r',ki) 
-B/>)N« (r,*0N« (r', fel ) 



(A3) 



(r,r'>R), 

G< 22 )(r,r» 

ik 2 



+K 



^ <H Kl + l) (l + m)\ 

1=1 m=0 \ ' ■> \ ' I 

x[c^(u;)M :im (r,k 2 )M :im (r\k 2 ) 
+ C?(Lj)N llm (r,k 2 )N :im (r',k 2 ) 
(r, r' < R), where 



UJ 



ki = -, k 2 = y/e(uj) — 
c c 



(A4) 



(A5) 



M and N represent TE and TM waves, respectively, 

777 ( sin \ 

M. Im (r,A0 = +-L-j l (kr)Pr (cose) i (m0)e e 
smO \cos/ 

. „ ,dPr(cos9) /cosV ,, . 

-^^^Hsin^' (A6) 



N = ; m (r, k) 



1(1 + 1] 



kr 

1 d[rj[(kr)] 



jt(kr)Pr (cos 0) Q (m(t>)e r 
dP] 71 (cos 6>) /cos 



kr dr 



dO \sin 



sm9 Vcos/ 



(m<f>)eg 



(A7) 



with ji(x) and P[ n (x) being respectively the spherical 
Bessel function of the first kind and the associated Leg- 
endre function. The superscript (1) in Eqs. flA2j) and 
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( [A3| ) indicates that in Eqs. ( |A6| ) and (A7) the spherical 
Bessel function ji (x) has to be replaced by the first-type 
spherical Hankel function /i|^(x). The coefficients Bf I,N 



and C ; M,Ar in Eqs. (A3) and ( |A3| ) are defined by 
[fgjrCgg)] ji(zi) - [gijfC^Q] J; (2:2) 



(A8) 



fHjj (^2) [fijj (zi )] ' - jj (zi) [z 2 ji {Z2j\ ' 
s(w)j l {z 3 )[z 1 h\ 1 \z 1 )]'-[z a ji(z2)]'hj 1 \z 1 ) ' 



[« 2 i I (^)]W 1) (*i)-i«(^)[^ft I (1) (^)]' 



(A10) 



; H^ (1) (z 2 )[z 1 / l | 1) (z 1 )]'-/ l | 1) (z 1 )[z 2 / l | 1) (z 2 )]' 



(w)ii(2 2 ) («0] ' - [wife)] (zi) 



where 

Zi — ki It 

Note that the relations 



(AH) 
(A12) 



1 fZ + nV 

dx P? (x)P r l (x) = S^TXTT^ ^ A13 ) 
.! (/-n)!(i + l/2) 



and 



/i, (1) (z) -> z-^xp^z-Ztt^-tt/^] if |z|^oo 



(A14) 



are valid |26|. 



APPENDIX B: SG FIELD RESONANCES 

For large value off = Z + l/2 the following asymptotic 
expansions are valid pa : 



J u (nx) 



exp[i/(tanha — a) 
V27wtanha 



u/c (coth a) 



fc=i 



(Bl) 



exp tanh/3)] 
^77^ tanh /3 



fe M fc (coth/3) 



fc=i 



(B2) 



J'Anx) 



sinh 2a 

47W 



exp[i/(tanh a — a)\ 



E 



Wfc(cotha) 



(B3) 



sinh 2/3 



y(/3-tanh/3) 



fc Wfe(cOth_a) 



fc=l 



[y^(x) - Neumann function], where 



CO 

x = —R, 

c 

V V 

cosh a = — , cosh/3 = — 
nx x 



(B4) 



(B5) 
(B6) 



and Uk and are given in Ref . [ p6[ | . To find the leading 
terms in Eqs. (0) and we thus may write 



J v (nx) 



I sinh 



«|, 



and 



H'„(x) Yl(x) 



H v {x) Y v {x) 



-| sinh /3 1 



to obtain, 



\/ v 2 — x 2 + \/ v 2 — ex 2 = 
for TE waves, and 



+ \/ v 2 — ex 2 = 



(B7) 



(B8) 



(B9) 



(BIO) 



for TM waves. Here we have used relationships flB6|), and 
we have assumed that x scales as v to discard the last 



term in Eq. (11) 



Obviously, Eq. ( B9 ) for TE waves has no solu tion, ex- 



cept for the trivial case of e= 1. Equation (BIO) for TM 
waves can be rewritten as 



(Bll) 



which just implies condition (|22|). The higher-order cor- 
rections can be obtained by writing 



x = vy 1 + 1 



7 lk 



k=l 



(B12) 



expanding all the quantities in Eq. ( |TT[ ) in powers of v 1 
and identifying the corresponding terms. 
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APPENDIX C: NEAR-SURFACE LIMIT 



Using the asymptotic Bessel-function expansion |2f 



Ju{z) 



YAz) 



oTT ) - (CI) 



>1), the coefficient B^{uj a ) 
in the form of 



Eq. (A9), can be given 



hf\k A rA) 



kATA 



c{u)a) - 1 



i(k A r A ) 3 (2l + l)e(u J A) 



1 



2/+1 



(C2) 



When the atom is located very close to the surface of the 
microsphere, i.e., r A >R, then from Eq. (C2) it follows 
that the series in Eq. (^) converges very slo wly. Hence, 
it is a good approximation to apply equation ( |C2| ) to the 
terms with small I as well. In this way, we derive 



1 c(u>a) - 



1=1 



h^\k A r A ) 

kATA 
1 1 



M e(u) A ) + 1 (Ar) 3 



(C3) 



Substitution of this expression into Eq. ( p8| ) then yields 
the leading term in Eq. ([4l|). The leading term in Eq. ( [42] ) 
can be derived in a similar fashion. 
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